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CHAPTER 1 


INTRODUCTION 

The phase-locked loop has long been recognized as a circuit with many impor- 
tant applications, and as such, the description of analog phase-locked loops (APLL's) 
has become well known as a large volume of material has been published to facilitate 
their use. In a survey of phase-locked loop development Gupta [ 11 ] has listed thir- 
teen books and over 120 papers which discuss APLL design and applications. In recent 
years, there has been an increasing use of various types of loops employing discrete 
elements. Among these have been hybrid loops that utilize both analog and digital 
circuitry. Newer loop realizations have been circuits composed entirely of digital 
elements, the digital phase-locked loop (DPLL). The importance of these types of 
configurations lies in the relative ease of design and construction and, of equal im- 
portance, the ease in which such circuits can be maintained. 

Unfortunately, the very attributes that make DPLL's attractive from a design 
and construction standpoint also contribute to difficulties in the theoretical analysis 
of DPLL operation and in this area the available literature is relatively thin. The 
first description of a general DPLL model was given by Reddy and Gupta [ 3 ] and 
further discussed by Gill and Gupta [ 10 ]. A more restricted but useful DPLL model 
for second order DPLL's has been presented by Weinburg and Liu [4 3. However for 
both of these general DPLL models, little has been done to describe the response of 
the DPLL underspecified ideal inputs. 

In the area of fading input analysis of DPLL's the available literature may be 
separated into two approaches. The first approach as employed by Holmes and 






Tegenelia [ 13 ], Weinburg and Liu [4 ], and Lee, Harington, and Cox [12 ], has 
been to model the DPLL configuration under study to operate linearly under some set 
of linearizing assumptions and then use classical techniques to determine the char- 
acteristics of the loop response. The second approach is more limited in that spe- 
cific loop configurations are assumed and are then analyzed by random walk tech- 
niques. This method was first used by Holmes C6] and has the advantage that no 
linearizing assumptions are required. This approach has also been used by Cessna 
and Levy [5], Yamamoto and Mori [14], and Ransom and Gupta [15] for other spe- 
cific loop configurations all of which employed resetting loop filters. 

In the following chapters, both the time response and the steady state noise 
response of first- and second-order DPLL's will be discussed in detail. The loop 
configuration used throughout will be based on the model described by Weinburg 
and Liu [4], hov/ever the results will not be restricted by the linearizing assumptions 
made in reference [4], Since the initial impetus for the study of the digital phase 
lock problem was provided by an application to an Omega navigation receiver, the 
practical realization of the DPLL"s considered will be presented and all examples 
and parameter selection will be based on this application. However, neither the 
analysis techniques developed to describe the DPLL°s operation nor the general 
hardware design presented are limited in any way to this application solely. 

In Chapter II, the APLL is described to provide a background for the phase lock prob- 
lem. In Chapter III, the response of the DPLL to specified ideal input signals is 
considered. In particular, new expressions are derived for frequency lock range 
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for firsf and second order DPLL’s and a new partial solution for the difference equa- 
tion describing second order DPLL operation is obtained. Chapters IV and V are de- 
voted to the steady state noise analysis of first and second order DPLL's respectively. 

In both cases, the loops are modeled as first order Markov chains and from this 
model the steady state phase error and mean time to phase lock are determined for 

an input signal corrupted by white gaussian noise. For the first order DPLL, a loop 

employing a non-resetting loop filter is analyzed and found to have essentially 

equivilant steady state phase error os for the loops using Random-Walk loop filters 

described in references [5], [14], and [15], but with less circuit complexity. 

Chapter V presents a new analysis for a second order DPLL, using the Markov chain 

model, that is valid for all values of input signal-to-noise ratio. The previous 

methods of analysis for second order DPLL's as described in references [4] and [13] 

have involved linearizing assumtions that have limited the usefulness of the analysis 

to signal-to-noise ratios greater than 0.0 dB. Experimental verification of the 

Markov chain model was desired to show validity of the noise analysis of Chapters 

IV and V and therefore, in Chapter VI, a description of the design and construction 

of a first- and second-order DPLL is given. The experimental test data used in 

Chapters IV and V were obtained from this loop design. 


CHAPTER II 


THE ANALOG PHASE- LOCKED LOOP 
A. Introducfion. The standard analog phase “locked loop as shown in 

Figure 2-1 has been widely described in the literature for a variety of applications. 
The following brief description of the APLL is presented to develop the background 
material necessary fora comparison with the DP LL operational characteristics to be 
developed in the following chapters. In particular, it 'vill be necessary to develop 
APLL response characteristics for both ideal and fading input signals. 

As can be seen from Figure 2-1 , the APLL consists only of a phase detector 
(multiplier), a linear filter, and a voltage-controlled oscillator (VCO) arranged 
in a feedback loop, in operation, the APLL is essentially a very narrowband filter 
whose characteristics are significantly dependent upon the type of linear filter em- 
ployed in the loop. The derivation that follows is similar to that by Viterbi [ 1] 
and Lindsey and Simon [2]. 


PHASE 

DETECTOR 
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B. The AFLL with Ideal Input . Consider the case of an ideal signal of 


the form. 


r(t) = y*2 A sin [ u t + 0(t) ] 

o 


(2-1) 


where A is the total power of the input signal and u is the quiescent frequency 


of the VCO. The input is angle modulated by 0(t) which is given by 


0(t) = Q t + 0 (t) 
o o 


( 2 - 2 ) 


where is defined as the frequency offset from the VCO quiescent frequency 
and 9^(t) is some function. The reference signal, y(t) at the VCO output can be 


expressed as. 


y(t) = cos + Ky^Q/e(T)dT] (2-3a) 

= -/2 K, cos[u t + ^(t)] (2-3b) 

I o 


2 

where is the total power of the reference signal, I*'yco'^ VCO 

with units rads/sec- volt, and 


"§‘(0 = Ky(~Q/e(T)dT 


(2-4) 


is the instanteous phase estimate of the reference signal. 

For the input and reference signals given, the error signal x(t) at the phase 
detector output is given by, 

x(t) = AK^ sin [ 0(t) -%) ] 

I m 
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(2-5) 


where the term involving 2o^t has been ignored since it will be removed by the 
loop filter/VCO combination and is the phase detector (multiplier) gain. If 
the initial conditions of the loop filter are zerO/ then the filter output e(t) can be 
expressed as, 

e(t) = A' K f(t - t) sin [0(t) -'§(t) ] dx (2-6) 

m o 

where f(x) is the impulse response of the linear filter. Defining the phase error 
(j)(t) for the loop as, 

<t>(0= 0(0 -e(t) (2-7) 

and substituting (2 -6) into (2-4) and taking the derivative gives, 

=-^ - AK f(t - t) sin <|)(t) dx (2-8) 

where K = ^m^VCO gain. Thus given some input phase 

function 0(t), the solution of the integro-differential equation of (2-8) for (j)(t) com- 
pletely describes the operation of the APLL. Equation (2-8) also suggests the stan- 
dard APLL model as given in Figure 2-2. 

C. The Linear APLL Model . The APLL model of Figure 2-2 can be 
simplified and more readily described mathematically to a good approximation 
whenever the phase error is less than 0.5 radians. In this case the approximation 

sin (})(t) ^ (j)(t) (2-9) 

can be made so that (2-8) becomes, 

^ - AK / f(t - x) ^ t ) dx (2-10) 

dt dt o 
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r 







Figure 2-2 . Standard APLL Model . 


Taking the Laplace transform of (2-11) with initial conditions assumed to be 
gives the frequency domain phase error as. 


s + AK F(s) 


0(s) 


(2 


Since <|)(s) = 9{s) “ 9(s), equation (2-11) leads directly to, 


A 

0(s) = 


AK F(s) 
s + AK F(s) 


0(s) 


^ H(s) e(s) 


where 


H(s) = 


AK F(s) 
s + AK F(s) 


(2 


zero 
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:-12a 

:-12b 

:-l3) 


is the closed loop transfer function of the linear APLL. Equation (2-1 2a) suggests 
the s-domain linear APLL model shown in Figure 2-3. 





0(s) 






AK F(s) 







1/i 





Figure 2-3. Frequency Domain Model of Linear APLL. 


Under the linearizing assumption of (2-9), the APLL response as described 
by (2-11) and (2-12a) can be readily analyzed by classical techniques and defi- 
nitions. Also, the APLL can be classified according to the characteristics of the 
open- loop transfer function G(s), 

F(s) 

G(s) = AK ^ — (2-14) 

In general for any control system, the order of the system is equal to the number of 
finite poles of the open- loop transfer funclion. Thus from (2-14) on APLL would 
be classified as a first-order system if there were no linear filter included in the 
loop. That is, F(s) = 1, Similarly, if, 

F(s) =1 + a/s (2-15) 

so that a signal-plus-integral loop filter is present in the APLL, then the APLL 
would be classified as a second-order system. 

Using classical techniques, consider the case of a first-order APLL with a 
frequency step input. This corresponds to, 
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I 


e(t) = Q t + e (t) 

o o 

which after toking the Laplace transform with 9^(t) a constant gives. 


0(s) = 


Q ,0 

0+0 


(2-17) 


Substituting (2-17) into C^-H) with F(s) = 1 gives 


(|)(s) = 


s + AK 




(2-18) 


From the final value theorem, the steady-state phase error for the APLL can be 
found as 


lim <^(t) = £^/ AK 
t -► 00 


(2-19) 


Thus the first-order APLL can achieve frequency synchronization but will not 
achieve perfect phase synchronization unless the frequency offeet between the VCO 
and the input signal is zero. 

Next consider the case in which the loop filter is the imperfect integrator of 
Figure 2-4. For this filter mechanization. 



Figure 2-4. Imperfect Integrator Loop Filter 


1 ,-^ . 


F(s) = 


1 + T^S 

1 + T^S 


where 


T, = (R, + R^) c 

T 2 = RjC 


( 2 - 20 ) 


( 2 - 21 ) 

( 2 - 22 ) 


Substituting (2-20) and (2-17) into C^-ll) gives. 


s (s + 1/t.) 

ct>(s) = ^ ! 

s + (AK T^/ + 1 /t^) s 



(2-23) 


and again applying the final value theorem gives the steady-state phase error as, 

a 

lim ({)(t) = (2-24) 

t “► 00 A 


which is equal to the stady-state phase error of the first-order loop attenuated by 
. Since in practice it is usually desired that the filter be of the signal-plus- 
integral type, will be a large number so that the pole of the loop filter is as 
near the origin as possible. Once again, this loop is capable of frequency syn- 
chronization but will track the phase with a constant offset. 

D. The APLL in Additive Noise. The effects of additive noise on APLL 
operation may be studied by considering an input signal of the form, 

r(t) = ^2 A sin [ u t + 0(t) ] + n.(t) (2-25) 

o I 

where n. (t) is the narrowband noise process given by 

n. (t) = [ X (t) cos u t - X (t) sin u t ] (2-26) 

I COSO 
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The terms x (t) and x (t) are assumed to be independent, stationary Gaussian 
c s 

white noise processes of zero mean and identical variances. Inclusion of the 
noise term gives a phase detector output of (referring to Figure 2-1) 


x(t) = AK^K^ sin (j<t) 

-K.K [X (t) cos <^t) - X (t) sin (})(t)] 
I m c s 


(2-27) 


where 


X (t) = X (t) cos 9(t) + X (t) sin 0(t) 
■j c s 


X (t) = X ft) sin 0(t) - X (t) cos 0(t) 
sc s 


(2-28a) 

(2-28b) 


Note that again the terms involving 2(o^t have been ignored. 

Following the same procedures as for the noiseless case, an i n teg ro- different! a I 
equation describing the APLL's operation is found to be, 

d(|)/dt = d0/dt - f(t - t) [ A sin <t>(T) (2-29) 

- X (t) cos (])(t) - X (t) sin d)(T) ] dx 
c s ^ 


As before, it is possible to represent the APLL by the simple model given in 
Figure 2-5. 

Applying the linearizing assumption, sln<jj~ <}) gives a simplified fading 
input model which can be analyzed by superposition. Thus it is only necessary 
that the noise terms be considered at this point since the effects of the determin- 
istic portion of the signal was considered earlier. Considering Just the noise 
term gives. 


-n- 




({.( 0 = ' 0(0 


(2-30) 


which will have a noise spectral density given as, 


S(j4)) - Sgdi)) 

KF(<o)/iu 

1 + 


2 

S (p) 
n 


(31a) 

(31b) 


where 3^t)) is the noise spectral density of the noise input. If the input noise is 

assumed to be white with single-sided density N then, 

o 


S<|)(u) - 


KF(^)/iu 

1 + 


N 


(2-32) 


Recalling the closed- loop transfer function of (2-13), the phase error spectral 
density can be rewritten as. 


N 

Sip(p) = j 

2A 



(2-33) 


The steady-state variance of the phase error then is easily found as, 

2 


N B, 
o L 




(2-34) 


where Bj^, the loop noise bandwidth, is defined as 


®L 2it 




00 


H(ju ) 


d(o 


(2-35) 


Thus once a loop filter has been specified, the loop noise bandwidth can be cal- 
culated and from this the operation of the APLL in the presence of noise is easily 
determined. 
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CHAPTER III 


THE DIGITAL PHASE-LOCKED LOOP 

A. Introduction . For the APLL of Section II.B, it was possible to derive 
an integro-differential equation in terms of the input signal phase and the phase 
error of the tracking loop, the solution of which completely describes the APLL's 
operation. However, unlike the APLL a standard form of DPLL has not as yet 
developed in the literature. This is due in part to the non-linear operation of 
virtually all DPLL components. Still it is possible, for a specific DPLL configura- 
tion, to develop a difference equation that describes the DPLL's operation. This 
difference equation for the DPLL is analogous to the integro-differential equation 
of the APLL. 

Difference equations have been derived for first and second-order 
DPLL's of a somewhat general configuration by Reddy and Gupta [3] and Wienburg 
and Liu [4]. However, the response characteristics of the DPLL's under considera- 
tion were largely investigated by iteration of the difference equations. This approach 
is basically a simulation of an ideal DPLL. 

In the following, since the DP LL under consideration here is a slightly 
modified version of those given in references [3] and [4], the describing differ- 
ence equations wilr be derived in detail. Further, it is very useful for a designer 
to have some knowledge of the response characteristics of a circuit without perform- 
ing a simulation. Therefore, through analysis of the describing difference equation, 
the partial DPLL response characteristics for specified inputs are derived. This 
includes determination of time required to achieve phase lock for some initial phase 
offset and frequency range for phase acquisition. 
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B. The Digital Phase- Locked Loop . The DPLL configuration under con- 
sideration here is shown in Figure 3-1 . For this loop the input is sampled at the 
positive-going zero crossing of the reference clock and the output of the phase 
detector is quantized to plus or minus one depending on the sign of the sample. In 
this manner, the phase difference measurement between the reference clock and the 
input signal is essentially reduced to a determination of phase lead or phase lag of 
the reference clock with respect to the input signal. The two forward paths repre- 
sent a filter in the sense that they alter the phase detector output to create a signal 
suitable for adjustment of the reference clock phase. Notice that the reference 
clock phase will take on a finite number of values determined by and A^. 








.. . - ;-,*r ,«»• 






Operation of the loop can best be described by considering the case of 

an ideal input signal of the form, 

s(t) = A sin [u t + 0(t)] (3-1) 

c o 

where 

0(t) = (u - u)t + e.(t) (3-2) 

o I 

defines a frequency offset plus a modulating phase function. The reference clock 
is given by, 

r(t) = sin [u t + 0 (t) 1 (3-3) 

o o 

which is assumed to be quantized to N distinct phase states so that (3-3) may be 
written as, 

r(t) - sin [u)^t + ^ (N - 2i + 1) 1 i = 1 , 2, * * *, N 

(3-4) 

For example if the phase state of the reference clock is 

i = N/4 (3-5) 

when the reference clock is given by 

r(t) = sin ( " 2 "" ' 

For the sampling example depicted by Figure 3-2, the input is sampled at 
the positive-going zero crossing of the reference clock giving a positive value for 
the first sample. Since this first sample is positive in sign, the output of the phase 
detector, X (1) is +1 . This phase detector output is then modified by the two for- 
ward paths to produce the signal, 

y(l) = “ - A 2 radians (3-7) 
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Input 



Figure 3-2. Waveform Sampling. 

which is the phase increment by which the reference clock is altered. Obviously, 
the values of and are related to the quantization of the loop. If at the 
second sample of the input signal (the second positive-going zero crossing of the 
reference clock) the reference clock still lags the input, then the reference clock 
will be altered in phase by, 

y(2) = “ A| - 2 A 2 radians (3-8) 

Similarly, if the lag condition still exists at the third sample, then the reference 
clock will be altered in phase by, 

y(3) = ~ A^ - 3 A 2 radians (3-9) 

This condition will continue until the j-th sample, at which time the reference 
clock phase leads the input signal phase and the reference clock will be updated 
in phase by, 
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(3-1 Oa) 


y(i) = + A, - A Z X(0 
' ^ i=l 

= + A^ - (i-DAj (3-lOb) 

Similarly/ if at the (j+1) -th sample the reference clock still leads the input then, 

y(j + 1) = + A^ + 2 A 2 - A 2 (i - 1) (3-11) 


The loop updating will continue in this manner until the DPLL achieves a lock 
condition which is characterized by 


y(0< 


A, + Aj 


(u - u) T(i) 


(3-12) 
.th 


for all successive samples with T(i) defined as the time interval between the i 
and (i - 1 )-th sample. 

For the analog phase-locked loop described In Section !!.B, it was 
possible to study loop characteristics from the solution of an integro-differential 
equation that modeled the loop's behavior. In an analogous manner, it is possible 
to develop a difference equation for the digital phase-locked loop that character- 
izes the loop phase error on a sample-to-sample basis. Because the loop phase 
error is determined in terms of loop sample number, it is also necessary to develop 
an equation describing the time of occurrence of a particular loop sample. The 
combined use of these two discrete equations then will completely model the opera- 
tion of the DPLL for an arbitrary phase input. The remainder of this section then 
will develop in detail these two describing equations. 

Proceeding in a manner similar to [3], the total phase of the input 
signal given by (3-1) is, 


P. (t) = + 0(f) 

I o 


(3-13) 
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and the total phase of the reference clock given by (3-3) i'S, 

S (t) = (0 t + 0 (t) (3-14) 

O o o 

However, the input signal can be sampled only at discrete times and the reference 
clock phase can exist only in specified states so that (3-13) and (3-14) are more 
properly written as, 

P.(k) = co^t(k) + 0(k) (3-15) 

and 

P (k) = u) t(k) + 0 (k) (3-16) 

o o o 

where for simplicity f [ t(k) ] is written as f(k) with k representing the sample 
number of the loop. 

The output of the phase detector is a sequence of values ±1 , the sign 
being determined by whether the phase of the reference clock leads or lags the 
phase of the input signal. Thus the phase detector output sequence, X(k), may be 
written as, 

X(k) = sgn [sin[p.(k)-p (k)] } (3-17a) 

I o 

= sgn {sin [ <)^k) ] } (3-1 7b) 

where 

#) = 0(k) - 0 (k) (3-18) 

o 

is defined as the phase error for the DPLL. Further, the phase of the reference 
clock is altered by a value Y(k) radians at the k-th sample so that the phase of the 
reference clock is given as, 

k-1 

0^(k) = I Y(i) + 0 q(o) (3-19) 
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where 9(o) is the initial phase of the reference clock which will be assumed to be 

zero. Substituting (3-2) and (3-19) into (3-18) gives the loop phase error as, 

k-1 

#) = (u -u)t(k) + 0.(k)-I Y(i) (3-20) 

' i=l 

As was stated previously, it is desired to develop a difference equation 
that is independent of time. Therefore, it is necessary to determine an expression 
for t(k) in (3-20) that is dependent only on the sample number k. First, define the 
value T (k) as the time interval between samples so that, 

T(k) = t(k)-t(k-l) (3-21) 

However, this can also be expressed as, 

T(k)=T-iYM) (3-22) 

ZTT 


2"ir 

where T= is the period of the quiescent reference clock. The time of occur- 

o 

rence of the k-th sample then is Just the sum of all T(i), i < k, or, 
k 

t(k) =I T(k) (3-23a) 

i=l 


T 

= kT-^ I Y(i) 


2tt . 


(3-2 3b) 


=1 


Substituting (3-23b) into (3-20) gives an expression for phase error, 
after algebraic manipulation, as, 

(i) k-1 

~ I 

O 1=1 


<{)(k) = 0.(k) + — 5 2irk I Y (i) 

I o_ 


(3-24) 


and the sample-to-samp!e difference in phase error is. 


# + 1)-#) = 0.(k + l)-9.(k)+2ir -— Y(k) 

I I u u 
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(3-25) 




' r ■ ^ 'V'. • ■■ 


The value Y(k) is defined as fhe phase increment of the reference clock for the 
k-th sample and will take on values determined by the gain constants and A 2 


Y(k) = A,X(k) + A, I X(i) 

^ i=] 


(3-26a) 


A^ sgn [ sin <|)(k)] + A 2 I sgn [sin (}>(i)] 

i=l 


(3-26b) 


Substituting (3-26b) Into (3-25) gives the phase error difference equation for the 
DPLL as. 


u “ tJ 

<b(k + 1 ) - <j)(k) = 0. (k + 1 ) - 9. (k) +2ir — 

1 I u 

o 


(3-27) 


{A. sgn [sin t(k)] + A- I sgn [ sin (j>({) ] } 

“o i=l 

To be complete, it is necessary to express (3-23b) in terms of the phase 
error by substituting (3-26b) into (3-23b) giving. 


t(k) = ^ [ sin (j)(i)] 

i=l 


(3-28) 


+ Ap r sqn [sin (jj(i)V, 


The DPLL then is completely described by the equation pair, (3-27) and (3-28), 
the former describing the phase error of the loop and the latter describing the 
discrete time ut which the phase of the reference clock may change state. 

C. First Order DPLL Time Response. A first order DPLL results for the 
system of Figure 3-1 when the gain term A 2 is set to zero. For this case (3-27) 
and (3-28) become, respectively. 
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(3-29) 


# + 1) - <^k) ■- 0.(k + 1) - e.(k) + 2u 


u - U) 


■ sgn 4<k) 


and 


k-1 


f(k) ^ ^ 50^9(0 

i= 1 

Note that for first order DPLL response the phase error is constrained to 

^(k) 


(3-30) 


V IT 


(3-31) 


so 


that 


sgn [sin <t(k)] - sgn (t(k) (3-32) 

As defined previously, the total number of phase states of the reference clock is N. 
Since at each sample of the input signal the reference clock phase is incremented 


by a value t A ^ , then 
N - 2ir 


so that (3-30) can be written as, 

T 

Kk) - kT - -jq- r sgn i|<i) 

hi 


(3-33) 


(3-34) 


Consider first the response of the first order DPLL to an input signal 
whose frequency is identical to that of the DPL L's reference clock but diffei’s 
initially in phase by some constant Pj. For this case, u ~ u and the input phase 
will be a constant with 9.(k ^ 1) P.(k) for all k which without loss of generality 

can be assumed to be zero. Thus, (3-29) becomes 

<j^ r 1 ) - d^k) “ ] ^9^ “^k) (3-35) 

with the boundary condition d(o) 0 j. For this simple case, the loop will be up- 
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dated by a value radians at each sample until the phase error changes sign 
at which point the quantized phase error will oscillate about the true phase of the 
input. An example response is given in Figure 3-3, It is evident that a lock con- 
dition occurs whenever 

<|)(k) I < A^ for all k > K (3-36) 

where K is given as , 

K = Integer [ 9j /A^ ] - 1 (3-37) 

The time of occurrence of the lock condition can then be found by evaluating (3-34) 
for K given by (3-37). Note that for all k < K the value of sgn (|j(k) will not change. 



Figure 3-3. Response of First-Order DPLL to Phase Step Input. 
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Next consider the case where the input signal is offset by a constant 
frequency from the first order DPLL reference clock. For this condition the phase 
error difference equation of (3-27) is given as, 

CJ - (i) 

(j<k + 1) - #) = 2 tt - — A, sgn 4<k) (3-38) 

o o 

where it has also been assumed that 0.(k) is equal to zero. In reference [ 3'^ an 
expression was obtained for frequency lock range by making the approximation 
<l}(k)~ 0 in fhe steady state which ignores the quantization of the loop. However, 
it is possible to include the effects of quantization in determining frequency lock 
range. Recalling the condition for first order DPLL. lock as given by (3-36), then 
an equivalent condition for lock is, 

I # + 1 ) - ct<k) I < 2 A ^ (3-39) 


which after some algebraic manipulation gives 

N-2 ^ (j ^ N+2 

N-1 - (. - N + ] 

o 


(3-40) 


This expression is plotted in Figure 3-4 and gives the normalized frequency range 
over which the first order DPLL can remain locked as a function of the number of 
states of the DPLL reference clock. Note that the undefined lower limit for N = 1 
is of no consequence since this would correspond to a DPLL with only one phase 
state. 


- 24 - 


1.5t 



Figure 3-4. Firsf Order DPLL Lock Bounds for Frequency Offset. 

D. Second Order DPLL Time Response. The time response for a first order 
DPLL has been characterized in the previous section for specified inputs. Notice 
however that the characterization of the time response did not involve the solution 
of the phase error difference equation because of the nonlinear nature of the equa- 
tion. Likewise, it is not possible to find a direct solution for the phase error differ- 
ence equation for a second order DPLL. Furthermore, no significant discussion of 
second order DPLL response has been found in the literature for the generalized 
model of Figure 3-1. The following, while not giving a complete solution, fills 
some of this void by determining valid and pertinent characteristics of the second 
order DPLL time response for specified input conditions. 
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As with the first order DPLL, the first case to be considered for the 
second order loop is an input signal with a constant initial phase offset 0| and 
zero frequency difference with respect to the reference clock. Also, the input 
signal phase will again be assumed to be zero radians and the constraint of (3-31) 
is applicable so that (3-27) becomes 

k 

(({k + 1) - 4(k) =-A, sgn ()3(k) - I sgn (|)(i) (3-41) 

i=l 

with <^ 0 ) = 9j. Since the total number of phase states of the reference clock is 
N , it is necessary that some relationship exist between A-j, A 2 / and N. For 

the second order loop it is assumed that A^ and A 2 are integrally related as, 

A^ = n A 2 (3-42) 

and the total number of phase states will be given as, 

N ---- 2tt/A^ - 2Tr/n A 2 (3-43) 

The typical phase error response versus the loop sample number as given 
by (3-41) is shown in Figure 3-5. While a complete solution of (3-41) is not per- 
formed it would be of interest to determine the sample number of each zero crossing 
of the phase error function and both the sample number and the value of each peak 
overshoot of the phase error function up to the occurrence of phase lock. A 
knowledge of these values not only characterizes pertinent parts of the phase error 
function but can be used in conjunction with one another to give the total time re- 
quired to achieve phase lock for an initial phase offset. Therefore, it will be 
necessary to determine the values k^ , k 2 , • . • , k.,kppkp 2 , • • • , kp., and 
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Clock Phase 



Figure 3-5 . Response of Second-Order DPLL to Phase Step Input. 


■'’osl' '•’os2' • • •1’osi 




OSI < 


A, + Aj 


(3-44) 


at which time a lock condition has been achieved at k.. 

I 

If the initial phase offset is assumed to be 0| > 0, then the first zero 
crossing of the phase error function occurs when 

#,) < 0 (3-45) 

This value of k^ may be found by considering the phase error for values of k < k^ 

which is given as , 


- -1 


#) = 0 - k A. - I i 

i=l 
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(3-46) 




(3-47) 


The first zero crossing of the phase error function occurs at 


kA^ + Aj I I > e, 

i=l 


However, 


k 

I i 
i=l 


k(k + 1 ) 


(3-48) 


which may be substituted into (3-47) to give (after rearranging terms), 

2 


2 2 A^ + A 2 

k + : k 


A, 


0 . > 0 


(3-49) 


The positive root solution for (3-49) evaluated for equality to zero gives the number 


of samples required for the first zero crossing of the phase error function as, 

0 . 


- (2 A, + A J 

k = ! + 

1 2 Ao 


2 A^ + A 2 ^ 2 


2 A, 


(3-50) 


After the k - th sample the sign of the phase detector output changes giving 
the phase error as , 

+ 1) = A, - A 2 (k, - 1) (3-51) 


which can be generalized to 


()<k^ + i) = (Kk^) + i A^ - — i (2k^ - i - 1) 


(3-52) 


for k, + i < k .. Defining k , as the sample for which the peak overshoot 

1 - pi pi 

occurs following the k. zero crossing, then it is obvious that, 

9(l<pi) < 9(kp; ^ P (3_53) 

Applying the constraint of (3-53) to (3-52) gives the sample at which peak over- 
shoot occurs as , 
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k . = 2k. - A,/A« - 1 
pi I 12 


= 2k^ - n - 1 


(3-54) 

(3-55) 


The vaiueof the peak overshoot is found by evaluating (3-52) at i=k ^ with 

k , determined from (3-55) giving, 

P' 

4<kp,) = <j<kp + (k^ - n - 1) r A^ - (k, + n) ) (3-56) 

It is important to note that at the peak overshoot point that the value contained 

in the summation block of Figure 3-1 will be zero. 

The partial response of the loop giving the zero crossings and the peak 

overshoot values for the phase error function can be found by repeated application 

of (3-50) and (3-56), Notice that since the summation term of the phase error 

function is zero at sample k , then k. , can be found from (3-50) where the 

pi 1+1 ' ' 

initial phase offset is given by Further, once the values of A^ and A 2 

have been specified it is a simple manner to determine the required response values 

2 

graphically. As an example, consider the case for which A ^ = 6 2 = — • 

For this case (3-50) reduces to. 


k^ =-1.5 +- (2.25 + 10.2 0|)‘ (3-57) 

and the phase error of the first zero crossing is determined from (3-46) as. 


(l<k) = 0| - 0,196 (k^ +- 3k^ ) 


(3-58) 


Equation (3-58) is in slope-intercept form with the intercept determined by (3-57). 

Notice however that the value of k determined by (3-57) applies only over a 

1 

certain range of 0, and therefore (3-58) is likewise applicable only over the same 
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range of 0|. Thus, when (3-58) is plofted in Figure 3-6 it takes on a saw-tooth 
characteristic. Substituting the loop gain values into (3-56) gives the value of the 
first peak overshoot as, 

0>(l<pl) = 4(l<i) + 0.098 (k^ - 2) (1 -kp) (3-59) 

which is also in slope-intercept form and is plotted in Figure 3-7. Notice that when . 
the zero crossing occurs at the second loop sample, the loop will be in phase lock. 

It is noted that the phase error at the occurrence of the first zero cross- 
ing of the phase error is plotted on the ordinate of Figure 3-6 and on the abscissa of 
Figure 3-7. Further, since the i-th zero crossing parameters are found by using the 
above procedure with the initial phase offset equal to the value of the (i-l)-th peak 
overshoot, then the abscissa variable of Figure 3-6 is equal to the ordinate variable of 
Figure 3-7. Therefore ifis possible to combine the results of these two plots. This is 
done in Figure 3-8 where (3-58) for 9| > 0 is plotted in the first quadrant and for 
0| < 0 in the third quadrant„ Similarly, (3-59) for a negative phase error at the 
zero crossing is plotted in the second quadrant and for a positive phase error at the 
zero crossing in the fourth quadrant. 

Using Figure 3-8 it is possible to determine the phase error character- 
istics during loop phase acquisition as follows. Assume that the initial phase offset 
of the loop is 0| = 2.75 radians as shown in Figure 3-8. The first zero crossing of 
the phase error function will occur at the fifth sample and the phase error will be 
-0.98 radians. The value of the first peak overshoot is the projection on the -0| 

axis of the intersection of <^k) = 0.98 radians with the k = 5 line in the second qua- 
drant which is equal to -2.04 radians. This value of first peak overshoot is then used 
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-0(k ) (Rad) 





fo find thaf flie second zero crossing occurs four samples offer fhe peak overshoof 
and has value 0.72 radians. The value of the second peak overshoof is seen fo be 
1 .3 radians. The value of fhe second peak overshoof is then used fo find fhe third 
zero crossing and fhe third peak overshoof value. This process is just continued in a 
clockwise manner until only two samples are required for fhe occurrence of fhe i-fh 
zero crossing at which time phase lock has occurred. Thus the phase function char- 
acteristics for any initial phase offset can easily be determined by Figure 3-8. 

Consider next the case where the input signal differs by a constant 
frequency from the quiescent frequency of the second order DPLL with 0.(k) = 0. 
Since for this case the condition of (3-31) is not assured during loop phase acquisi- 
t’on the phase error difference equation is given by 

(f(k + 1) - tj)(k) = 2 )t ( 1 - u/u ) - ^ ~ — sgn [ sin 0(k)] (3-60) 

o 1 to 

k ° 

- u/(j T sgn [sin0(k)] 

° i= 1 

with (j)(0) = 0|. The inclusion of the constant term for the frequency difference and 
the sine function complicates (3-60) to the point that even a partial solution as was 
performed for the case of zero frequency difference is not possible. However, it is 
possible to study the stability of the second order DPLL by means of an incremental 
phase plane portrait. In the incremental phase plane, <jT(k + 1) - <h(k) is plotted as 
the ordinate and (^k + 1) is plotted as the abscissa while k is a variable parameter. 
Then for any k, f <j3(k + 1), 4>(k + 1) - 4(k)] describes the state of the system and a 
stable region of operation is found when, for any given initial conditions, 

r# + 1), (Kk + 1) ~(j5(k)] r2iTi, 0.0] i = 0 , ± 1 , x2 . . . (3-61) 


as 


k -*-co 
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For the case of = A 2 = ir/32 and 0| = ir/4 Incremental 
phase plane portraits have been plotted in Figures 3-9 and 3-10 for various values of 
(ii/ From Figure 3-9 it is seen that for < 1/ but sufficiently large, that it 

is possible to achieve phase lock without cycle slipping. However, as decreases/ 

the loop will achieve phase lock after slipping an increasing number of cycles. Notice 
that as the number of cycles slipped increases the length of time required to achieve 
phase lock also increases. Similarly, from Figure 3-10 it is seen that for > 1/ 

but sufficiently small, that it is possible to achieve phase lock without cycle slipping 
but as (o/ to increases the number of cycles slipped prior to phase 'ock increases. 



J- - it/2 

Figure 3-9. Incremental Phase Plane Portrait forO, = Tr/4andto/ to < 1. 
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Figure 3-10. Incremeni-al Phase Plane Porfraif for 0 - n^4 and u/u 


For some applicafions it is necessary that phase lock be achieved in 

minimum time so that it is necessary that no cycle slipping occur during phase acqui 

sition. Consider first the case of 0, >0 and u /cj < 1 , so that the phase error 

I o 

can be expressed as 


&{k + 1 ) - 0, + 2tt ( 1 - u .''u ) k - Q /u A , k 
1 o' o 1 


(3-62' 


/ . k(k + 1) 

- u/o Ao 

o 2 J 

up to either the first zero crossing of the phase error function or the first cycle slip 
s nee sgn f sin <b(k) ] - 1 until either of these occurrences. Equation (3-62) can 


be rewritten as 


(b(k +1) = 0j - A.J k-A 


k(k - 1) 


(3-63; 


V 'here 


A, ~ -2tt+ l)/(j (27t4 a,) 

1 o 1 


(3-64; 


A« - (j/ci) Ao 
z o z 


(3-65; 
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iinH for a cycle slip to occur 
(^k * 1) TT 
From (3-63) and (3-66) 


r \ f \ 

, L_^ k (0 -n) 

A 2 A 2 


(3-66) 


(3-67) 


Since k Is the sample number of the occurrence of a cycle slip during phase acquisi- 
tion, k must be a positive, real value. Since k must be real then 


2 A^ * A 2 

2 a' 


2 (0 it) > 0 


(3-68) 


and, since 0, it, fork to be positive then. 


2 A ^ 1 A 2 


2 A „ 


(3-69) 


f'or (3-68), the worst case for cycle slip occurs as 0, -► it, therefore if 0, ^ it - 8 


and A 


TT ”32 the constraint of (3-68) gives 


u 'u 0.971 for 8 0.01 

o - 


(3-70) 


while the constraint of (3-69) gives 

o _ 0.9:77 

Therefore, k will be a positive real value for indicating a cycle slip will occur for 

all values of u cj •. 0.971 , A similar procedure can be followed for 

- TT •' 8 ' 0 and u o 1 which for the same case of A,= A« ~ ir/32 

- I - o 12 

will give a c.ycle slip for 


- 36 - 


(3-72) 


(j / (j 1 . 03 , 5 - 0 . 0 1 

o - 

Thus once d ^ and A2 are defined, If is possible for the 2nd order DPLL to define 
the frequency offset range over which phase acquisition can occur without cycle 
slipping. 


CHAPTER IV 


FIRST ORDER DPLL WITH ADDITIVE NOISE INPUT 
A. IntroducHon . In f-he preceeding chapter partial response characteristics 

for a first- and second-order DPLL with ideal input were developed. As for the case 
of an ideal input, little has appeared In the literature concerning the characteristics 
of DPLL operation with a noisy input. For the most part any analysis that has appeared 
in the literatuie for stochastic inputs has involved linearizing assumptions on the DPLL's 
operation. Two notable exceptions have been the DPLL configurations studied by 
Cessna and Levy f 5] and Holmes [ 6 ] where random walk techniques were used to 
d(itermine the statistical characteristics of the loop. For the DPLL of reference fS 
the first order loop utilized resettable low pass digital filters in the phase correction 
path while the DPLL in reference [ 6 ] did not utilize a filter. A DPLL configuration 
that uses a non-resetting filter is analyzed in the following without linearizing assump- 
tions using a Markov chain model that achieves the same statistical performance of that 
described in [ 5 1 but with less physical complexity. 


B. First-Order DPLL Configuration. The specific loop configuration consid- 
ered here is shown in its implementation form in Figure 4-1 and is a slightly modified 
version of the generalized DPLL model shown in conceptual form in Figure 3-1 with 
A 2 ■ 0- This loop configuration is that of the Ohio University MAPLL [ 7^ with 
the exception that the loop is assumed to operate continuously instead of in a gated 
manner. Assuming for the moment that M ‘ 1 in Figure 4-1 , then this loop implemen- 
tation is easily seen to be identical ro the model of Figure 3-1 with A 2 = 0 since at 

each sample of the input signal the phase of the reference clock will either be advanced 

2 

or retarded by A ^ radians depending on whether the reference leads or lags the 

input signal „ 
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An improvement can be made in the loop's operation for a fading input by 
the addition of the divide-by M up/down counter proceeding the divide-by N up/ 
down counter. It is obvious that to change the phase of the reference clock it is 
necessary for the divide-by M counter to cycle through its M distinct states to either 
an overflow of underflow condition. Thus, while the phase output takes on N distinct 
states, the loop itself has MxN distinct states. Further, the divide-by M counter 
constitutes a digital low-pass filter whose transfer function is found In Appendix D as, 

1 sin M — ^ „ -i (JTT 


“ M --7^ 

■-T/U 


sm uir/ 


e ■ — (M-1) 
u 

c 


(4-1) 


c 

It is seen from the plots of the magnitude-squared function of (4-1) in Appendix D 
that as M increases, the bandwidth of the filter decreases. 

A state diagram for the loop is given in Figure 4-2 where the values of p. 
and q. are the probabilities associated with the indicated state changes. Several 
things are worth noting about this state diagram. First, for any present state, when a 
new sample is taken a new state will result; and second, the new state will always be 
adjacent to the previous state. Also, for a given reference clock state the transitions 
occur uniformly with time but, following a transition from one reference clock state 
to another, the time interval to the next sample is either longer or shorter than the 
time interval between the previous two samples depending on whether the reference 
clock phase was advanced or retarded as it passed from the previous reference clock 
state to the present reference clock state. 

If the phase samples are independent and the non-uniform sampling interval is 
ignored, then the DPLL given by the state diagram of Figure 4-2 can be approximated 
by a first-order Markov chain, the properties of which are outlined in Appendix A. 
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Or, 


rp^ [a] = [a] (4-3) 

where [ p] is called the matrix of transition probabilities. The non-trivial solution 
to this system of homogeneous, linear equation is the one that satisfies the total proba- 
bility relationship , 


N M 

T T a.. = 1 
i=l '' 


(4-4) 


and, from the characteristics of Markov chains, represents the steady state probability 
for the loop states. That is, at any given instant of time the probability of observing 
the loop in state s.. is given by a., . 

The steady state probabilities for the reference clock states, s. are given 


by. 


M 

S. = r a.. i = l, 2, . . . N 

I . 1 U 
i=*i 


(4-5) 


Once the steady state reference clock state probabilities are known, it is possible to 


find the variance of the reference clock phase from 

N-1 


2 - 1 T 

‘’® - FT 


S. 


1=0 


i+1 


N 


( 1 - N + 2i) 


(4-6) 


where it has been assumed that the phase of the input signal is a constant 0. radians 
corrupted by noise. 

Another quantity useful for the evaluation of a DPLL with a fading input 
is the mean time to lock-up for some initial phase offset. This can be evaluated by 
considering the loop transient response in terms of the classical ruin problem (Feller 
[ 8]) for a random walk. The mean time to lock for an initial phase offset then would 
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be equivalent to a determination of the expected duration of the game in the classi- 


cal ruin problem. If T^. , is defined as the expected time required to reach a mini- 
mum phase error when the initial state of the loop is s,. then the T . . must satisfy 

l| 01,1 

the difference equation, 


T. . — 1+p. T, . i+q, T.. + 1 

oi , I I oi , I - I I oi ,1 

for J 1 , M and with the boundary conditions 


(4-7) 


Tok,r\.I,i=° i = l,2,...M (4-8) 

where reference clock states k and k +1 represent the minimum phase error states. 
For the input signal assumed to be 0. radians then the system of non-homogeneous , 
linear simultaneous aquations results. 


- 42 - 







c. 


Matrix of Transition Probabilities, 


Consider the case where the input to the DPLL is of the form; 
s (t) = s(t) -t- n(t) (4- 10a) 


= A cos (c t + i|)) + x(t) cos u t + y(t) sin u t 
c c Co 


(4- 10b) 


where A is the carrier amplitude and x(t) and y(t) are zero mean independent gaussian 
c 

distributed random processes of bandwidth B and variance <7^2 = a2 = cr2. That is, the 
input to the DPLL is some signal plus narrowband noise. The input can also be written 
in the form: 


wh : 


ire 


s (t) = x'(t) cosu t + y'(t)sinu t 
r c c 


x'(t) = x(t) + A cos ijj 


(4-11) 

(4- 12a) 
(4- 12b) 


y'(t) = y(t) + A^sin tjj 

The positive going zero crossing is always assumed to be the correct phase of the signal 
s(t). That is, the reference clock is always assumed to be in phase lock with the signal 
s(t), and the error signal generated by the phase detector is used to tell the loop differ- 
ently. Therefore, the loop always assumes the samples of the incoming signal occur at: 

(4-13) 


so 


that 


u t. = (1 + 2i)Tr/2 i = 0, 1, 2, . . . 
c I 


cos (j t = 0 
c 


giving 


sin u) t = 1 
c 


s^(t) = y(t) + A^ sin ijj 


(4- 14a) 
(4- 14b) 

(4-15) 
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where ijj is flie phase difference between s(t) and the reference clock. Note that as 
the loop approaches lock, i() approaches 0. Also, the reference clock can take on 


only N distinct values so that tp too can take on only N distinct values ip. , i = 1 , 
2, . . . , N. From earlier y(t) is a gaussian distributed random process so that the 
probability density function for s_^(t) is: 

1 


P(s ) = 


1 (y + A sin ip.)' 


/ 


2it I 


T 


(4-16) 


From the state diagram of Figure 3, the p. s are the probability that the sampled 


value of s (t) is less than zero and can be found from; 


p. = /° P (s I Ip.) ds 
r, ' r ^1 r 

-oo 


(4- 17a) 


Also, 


^2Tra^ 


y + A sin ip . 

_ i r ^ ' -I 

o ~2 [ j 

J e o 


GO 


dy 


= probability that s^(t) ^ 0 


(4- 17b) 


(4- 17c) 


q. = 1 - p. (4-18a) 

= probability that s^(t) > 0 (4- 18b) 

D. Application of the Markov Chain Model . The phase error and transient 
response of the DPLL given in Figure 4-1 can be determined from the solution of 
(4-2) and (4-9), respectively, once the values of M and N have been specified. To 
perform the solution to these two sets of equations, the three Fortran computer pro- 
grams given in Appendix B were written. The first program, PBSTGEN, is used to 

set up the matrix of transition probabilities for the systems of (4-2) and (4-9) once 
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the values of M and N are given. The second program given in Appendix B, 

PBDPL2, is used to solve for the steady state probabilities of the loop states and from 
this calculates the steady state probabilities of the reference clock states and finally 
the variance of the phase error for specified input signal-to-noise ratios. The final 
program given in Appendix B, PBDPL2T, is used to determine the solution to the 
transient response system of equations given in (4-9). Both PBDPL2 and PBDPL2T use 
a successive approximation technique to determine the solution to the system of equa- 
tions. Note, all three programs were written to be used both for the first-order DPLL 
considered in this chapter and also for the second-order DPLL which will be given in 
the next chapter. 

The steady state probabilities for the reference clock states were found for 
various values of M and signal-to-noise ratio with N=64 and are plotted in the terms 
of the phase error probability density function in Figures 4-3 and 4-4. The value of 
N = 64 wa chosen to be consistent with the application to an Omega navigation re- 
ceiver as ou'lined in [7]. Note that for all cases, the value of N and M will be 
assumed to be a power of two so that a practical realization of the counters can be 
achieved with a standard binary counter. In Figure 4-3, the envelope of the proba- 
bility density function for the phase error is shown for M = 1 and signal-to-noise 
ratios of 20. , 0.0, -20. , and 40. dB. Notice that as the signal-to-noise ratio de- 
creases, the envelope of the density function approaches that of a uniform phase error 
distribution as would be expected. In Figure 4-4, the envelope of the density func- 
tion is shown for a constant signal-to-noise ratio of -20, dB while M takes on values 
1,2,4, and 8. From this plot the effects of the digital low-pass filter (the divide- 
by M counter) can be seen since as the bandwidth of the filter decreases, the den- 
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sity functions tends to flatten. 


Once the steady state probabilities of the reference clock states are deter- 
mined the phase error variance can be easily calculated. In Figure 4-5, the standard 
deviation of the phase error is plotted for M = 1 , 2, 4, and 8 versus noise-to-signal 
ratio. For low values of noise-to-signal ratio these curves are asymptotic to a value 
fixed by N, the quantization level of the reference clock. For high noise-to-signal 
ratios, these curves are asymptotic to the standard deviation of a uniform phase dis- 
tribution. Note also that as the value of M increases by a power of two that there 
is an approximate 3.0 dB increase in the performance of the loop. 

The mean time to lock in terms of the number of samples required to 
achieve phase lock for some initial phase offset was determined for M = 1 , 2 , and 4 
for various values of signal-to-noise ratio. Note that phase lock following some 
initial phase offset is defined as the first occurrence of the minimum phase error con- 
dition as determined from the solution of (4-9). Since the phase of the reference 
clock is assumed to take on values from -tt to ir, a minimum phase error condition 
would be reached for M = 2 when the loop reached states s ^2 i / 5^2 2 ' ^33 1 ' 

S 33 2 for N = 64 and the input signal a constant 0.0 radians. However, for the 
case considered here, a more stringent lock condition was imposed in that only states 
s ^2 2 ®33 ] '^sre considered to constitute phase lock. For this case the bound- 
ary conditions for (4-9) become 2 ^o 33 1 zero. Similarly, for 

the case of M = 4, only the states s_„ ^ arid s„„ , were considered to constitute 

' 32,4 33,1 

phase lock. The solution to (4-8) under the above conditions are plotted in Figures 
4-6, 4-7, and 4-8. As expected, for a given initial phase offset, the mean number 
of samples to the first occurrence of a phase lock condition increases as the signal- 
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to-noise ratio decreases for all three plots. Notice that for a 20. dB signal-to-noise 
ratio that the expected time to lock is very nearly equal to the time required for lock 
under ideal conditions. A cross comparison between Figures 4 - 6 , 4 - 7 , and 4-8 indi- 
cates that the expected time to achieve phase lock increases as the value of M in- 
creases as would be required since more loop states must be traversed to achieve phase 
lock . 

In addition to the theoretical values, Figures 3-5 through 4-8 also include 
experimental data taken from a hardware realization of the first-order DPLL. This 
hardware realization will be discussed in Chapter VI. 
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Figure 4-1. First-Order DPLL Implementation. 

Note; i is Present Value of -s- N 
Up/Dn Counter. 
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V. SECOND ORDER DPLL WITH ADDITIVE NOISE INPUT 


A. Infroduction . Second order DPLL configurafions have been analyzed by 
Holmes and Tegenelia [13] and Weinburg and Liu [ 4] for an addifive whife gaus- 
.ian noise input. However, in the case of reference [ 13 ], the loop was modeled 
js a linear system and then analyzed by classical techniques. The results presented 
for this linear model showed close agreement with experimental data for signal-to- 
noise ratios greater than O.OdB. In reference [ 4] , the steady state phase error for 
the second order DPLL was determined from a solution of the Chapman-Kolmogrov 
equation in the z-domain under the assumption of small values of steady state phase 
error. Again, the usefulness of this solution is limited by the input signal-to-noise 
ratio. In this chapter, the author presents a method of analysis for a second order 
DPLL that is not limited by the input signai-to-noise ratio. 

In the previous chapter, the analysis of a first order DPLL was performed 
by modeling the loop as a first order Markov chain. In the following, the idea and 
methods used in the previous chapter will be modified to allow a similar approach to 
be used for the analysis of second order DPLL's. It will be shown that a second order 
DPLL can be modeled as a first order Markov chain with alternatives and that these 
alternatives themselves can be thought of as states in a first order Markov chain. 

The steady state distribution of the Markov chain alternatives can be determined and 
Tom this distribution it is possible to find the steady state phase error of the DPLL. 
The transient response of the loop is also determined in a similar manner. 

Be Second Order DPLL Configuration . The second order DPLL configuration 
considered in this paper is shown in its implementation form in Figure 5-1 . The loop 
is identical to that of Figure 4-1 with the exception of the addition of the divide-by 
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L up/down counter and the K-bit binary adder whereK = log„ (MxN). The in- 

k ^ 

elusion of the divide-by L counter provides the I X (I) function of the gen- 

^ I =1 

erallzed model of Figure 3-1. Note however that for a practical Implementation, 
the dIvIde-by L counter cannot be allowed to either overflow or underflow since 
this would have the effect of resetting the summation value to zero. Therefore, the 
divIde-by L counter Is structured so ttiat It will saturate at values of i L. 

Operation of the loop Is as follows. The Input signal Is sampled at the 
positive-going zero crossing of the reference clock to determine whether the refer- 
ence clock leads or lags the Input signal. If a phase lead Is detected, the dlvlde- 
by M and dIvIde-by L counters are Incremented by one and If phase lag Is detected, 
the dIvIde-by M and divide-by L counters are decremented by one. Following this, 
the value of the divIde-by M and dIvide-by N counter combination Is added to the 
value of the dIvIde-by L counter and the result loaded Into the dIvide-by M and 
dIvIde-by N counter combination. Finally, the phase of the reference clock Is up- 
dated to reflect the value contained In the dIvIde-by N counter after the load has 
occurred . 


Referring to Figure 3-1 , the value of A ^ will In this case still be given 
by A ^ = -jq- and the value of A 2 is determined by the bit In the dIvIde-by M and 
dIvIde-by N counter combination to which the feast significant bit of the dIvIde-by 
L counter is added. For example. If the binary adder is configured so that the least 
significant bit of the divide-by L counter adds to the least significant bit of the 
divIde-by N counter then = A 2 » However, If the least significant bit of the 
divIde-by L counter adds to the most significant bit of the divide-by M counter, 
then A^ =2 Note that this leads to a more limited relationship for A^ and A 2 
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than given in (3-42) since for this configuration, 

= 2* A2 i - 1, 2, . . . (5-1) 

Note also that a first order DPLL results from this configuration if the divide-by L 
counter is configured to saturate in its zero state. 

For the first order DPLL, the loop states were defined by the combination 
of the divide-by M and divide-by N counters giving a total of MxN loop states. 

For the second order DPLL, the loop states are defined by the value contained in 
the combination of the divide-by M and divide-by N counters following the loading 
of theK-bit binary adder output so that once again there are a total of MxN loop 
states. However, for the second order DPLL the loop state transitions are dependent 
not only on the present output of the phase detector but also upon the value con- 
tained in the divide-by L counter. Thus the value of the divide-by L counter can 
be thought of as providing alternative loop state to loop state transitions for a given 
phase detector output. The loop can then be modeled as a first order Markov chain 
with alternative state transition vectors. The characteristics of the Markov chain 
with alternatives is discussed in Appendix A. 

C. Markov Chain Model of the Second Order DPLL . Use of the Markov 
chain with alternatives for the modeling of the second order DPLL configuration con- 
sidered here will be shown by means of an example. For this, consider the DPLL of 
Figure 5-1 with M = 2, N = 4, and L = 3 so that there are eight loop states which 
define the states of a first order Markov chain. However the divide-by L counter 
can take on seven values so that for each loop state there are seven possible alterna- 
tive loop state transitions. A state diagram of the Markov chain model for the DPLL 

under consideration is shown in Figure 5-2. In this figure, the loop states are situ- 
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ated horizontally while the alternatives associated with each loop states are situ- 
ated perpendicularly Thus the top row labeled alternative 1 gives the state trans- 
itions for the divide-by L counter having a value 4-3. Similarly alternative 2 applies 
for the divide-by L. counter having value 2, alternative 3 is for the divide-by L 
counter having value 1 and so on. The loop states are plotted such that the first 
column is for a reference clock state of one with the value of the divide-by M 
counter equal to zero while the second column is for reference clock state one with 
the value of the divide-by M counter equal to one. Similarly, the third column is 
for reference clock state two with a value of zero contained in the divide-by M 
counter and so on . 

Since there are eight possible loop states with each loop state having 
seven possible alternative actions, there are a total of 56 alternative state vectors 
in the matrix of transition probabilities as given in Appendix A. However, for each 
loop state, the alternative vector to be used is uniquely defined by the value of the 
divide-by L counter so that the alternative vectors themselves can be thought of as 
states in a first order Markov chain. That is, instead of considering the transition 
from loop state to loop state, the transitions from loop alternative to loop alternative 
are considered. Thus in the state diagram of Figure 5-2, the loop alternatives are 
successively numbered s^, $ 2 / . . . and the possible transition from loop alterna- 
tive to loop alternative is assigned a probability as indicated by the directed arrows. 
The values of the indicated probabilities are dependent only upon the reference 
clock state associated with each loop alternative and the method of determining their 
value was discussed in Chapter IV. C. 
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To show l-hal- the state diagram of Figure 5-2 will achieve a phase lock 
condition, consider the case of an ideal input sigiva!. For this case the transition 
probabilities as given by (4-16b) and (4-17a) are, 

and 

P3 = P4 q] = q2 = 0 ( 5 - 3 ) 

If the loop is assumed to initially be in alternative state s^, then on successive 
samples of the input the path traced through the state diagram will be 

^4 - hz - ®37 ^38 ®32 - ®19 - ^25 - ®38 

so that when phase lock is achieved, the loop will continuously cycle through the 
closed path 

®38 ^32 ®19 ^^25 ^ ^38 ' ‘ ' 

Inspection of Figure 5-2 will shov/ that the same state cycle will be obtained for 
any initial alternative state. 

Once the state diagram for the alternative loop states is defined, it is 
clear that the steady state phase error for the loop may be obtained by the same 
method as was used for the first order DPLL in Chapter IV. For the example of 
Figure 5-2, if the steady state probability of occupancy of alternative state s. is 
denoted as a., then for state s, , 

I I 

“ P, =29 + q, O, ( 5 - 4 ) 

Similarly for alternative state 

“2 " Pi °29 + 'll “3 
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Obviously fhere exists a similar equation for each of the 56 alternative loop sta*-es 
so that Hiere exists a system of homogeneous, linear equations similar to (4-1) which 
can be written as. 


1 


Q 
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°2 
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°56 
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^56 


(5-6) 


where [A]^, ... is the stochastic matrix of transition probabilities for the alterna- 

56 X o6 

tive loop states. The nontrivial solution to (5-6) which satisfies the constraint, 


56 

r a. ] (5-7) 

. i I ' 

gives the steady state probabilities for the alternative loop states. As for the case of 
the first order DPLL, the steady state probabilities of the reference clock states, S., 

I 

con be determined from the steady state probabilities of the alternative loop states 


as. 
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I a. 
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I a. 


S 


3 


42 

I a. 

• O /-\ ^ 


(5-8c) 


(5-8b) 


(5-8c) 
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56 

I a. 

i 43 ' 


(5-8d) 


The variance of fhe steady state phase error con then be determined from the distribu- 
tion of S, by (4-5). This method of finding the steady state phase error of the loop can 
obviously be extended to any value of M, N, and L. 

The mean time to phase lock for some initial phase offset can also be deter- 
mined by the method of Chapter IV. If T^. is defined as the mean number of samples 
required to the first occurrence of phase lock for an initial alternative loop state s. 
then for state s ^ , 

^ol '^3 ^o29 P 3 ^o30 ^ 

A similar equation can be written for each of the alternative loop states with the 
exception of the alternative states which define phase lock for the loop. For the ex- 
ample of Figure 5-2, the alternative loop states defining phase lock were found earlier 
to be s.Q, s„ , s and s„ . For these states the loop is initially in phase lock so that, 

1 7 Z.O Jc 

^ol9 ^o25 ^o32 \)38 ’ 

Thus there exists a system of non -homogeneous , linear equations similar to (4-10) 
which can be written as. 
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where [T] is 56 x 56 matrix involving p. and q,. As with the steady state phase 
error, this method of determining the mean number of samples to phase lock can be 
extended for arbitrary values of M, N, and L. 

D. Application of the Markov Chain Model . The use of the Markov chain 
utilizing the alternative loop states can be generalized for arbitrary values of M, N, 
and L os has been done via the Fortran programs given in Appendix B. Use of these 
programs has been outlined in Chapter IV. C- For all coses considered it was assumed 
that and N =64. Further it was found in [9] that the phase error de- 

graded rapidly as L increased, so for all cases considered L was limited to values of 
0, 1 , or 3. 

Using PBDPL2 in Appendix B, the standard deviation of the steady state 
phase error was determined for L ~ 0, 1, and 3 for various signal-to noise ratios. 
Figures 5-3, 5-4, and 5-5 give the phase error for values of M = 1, 2, and 4, re- 
spectively. As can be seen from all three plots there is an approximate 20. db de- 
gradation in loop performance for L = 3 as opposed to a first order loop (S = 0). 

Using PBDPL2T given in Appendix B, the transient response of the second 
order loop was determined for M = 1,2, and 4, L = 1 and 3 and signal-to-noise 
ratios of -20., 0., and 20. db. The results are plotted in Figures 5-6 through 5-11 
where for all piots the value of the divide by L counter was zero at the initial phase 
offset. When the results are compared to those for the first order loop given in 
Figures 4-6, 4-7, and 4-8, it is seen that for identical values of M and signal-to- 
noise ratios, that the second order DPLL achieves lock In a lesser expected time for 
initial phase offsets greater than ir/8 radians. However for smaller initial phase off- 
sets, the first order DPLL actually achieves phase lock in a lesser expected time than 
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the second order DPLL. 
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Figure 5-6. Loop Transient Response 
M = l, N=64, L = 1 . 
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Figure 5-9. Loop Transient Response 
M = 2, N=64, L=3. 





CHAPTER VI 


FIRST AND SECOND ORDER DPLL DESIGN 

A. Infroducfion . To experimcnfally dernonsfrotc- fhe validity of the Markov 

chain model for the first and second order DPLL's described in Chapters IV and V, a 
hardware DPLL was designed, constructed and tested for additive noise inputs. The 
hardware loop was designed so thot either first or second order operation was selec- 
table. Also, a primary design objective for the hardware DPLL was to allow for easy 
alteration of important loop parameters to facilitate testing of the loop under various 
conditions. Therefore, the design presented is not intended to incorporate minimal 
components. The general range over which the parameters could be varied was limi- 
ted to values thought to be useful for an application to an Omega navigation receiver 
as were the parameters used it} Chapters IV and V. It should be emphasized however, 
that neither the theory nor the general DPLL design presented in this paper are in any 
way limited to this particular application. 

B. DFLL Dtsign and Testing. The block diagram tor the DPLL design is shown 
in Figure 6-1 . All register lengths shown in the block diagram indicate the maximum 
values and during testing the actual register lengths were altered to verify the effects 
of various loop parameters upon the DPLL's performance. The functional groupings of 
components in relation to the DPLL given in Figure 5-1 are indicated for comparison. 

Referring to Figure 6-1 , the phase detector will sample the binary quan- 
tized incoming signal at a frequency f and then output a count-down signal if the 
sample is a high level or a count-up signal if the sample is a low level . The count- 
up or count-down signal Is then applied to both the divide-by L saturating counter 
and the divide-by (MxN) counter. Inhibiting logic is included with the divide-by 

L counter so that the counter wi 1 1 saturate at sel ectable values of T (2 - 1 ) , 
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Figure 6-1 . Hardware DPLL Block Diagram 























i 1 , 2, . . , , 7. The sample command also initiates the control logic so that the 
new value of the divide-by fMxN) counter (following the count-up or count-down 
signal) is loaded into the 12-bit buffer by means of the LOAD2 signal. After set- 
tling, the output of the 12-bit adder will be the sum of the divide-by L saturating 
counter and the divide-by (MxN) counter. This value is then loaded into the 
divide-by (MxN) counter by means of the LOADl signal. The value of the N most 
significant bits of the divide-by (MxN) counter now represents the phase estimate 
of the DPLL. To establish the variable phase reference clock, this phase estimate is 
compared to the value of a divide-by N counter being clocked at a rate Nxf using 
a log 2 N - bit binary magnitude comparator. Upon coincidence of the two input words, 
the magnitude comparator output takes on value ONE. Note that this pulse output 
occurs at an f^ rate. Notice also that this loop will operate in the first order mode 
simply by inhibiting the LOADl signal to the divide-by (MxN) counter. Detailed 

schematics fcr the DPLL may be found in Appendix C. 

The primary objective for construction and testing of the hardware DPLL 

was the verification of the theoretical data obtained from the Markov chain model in 
Chapters IV and V. As such, it was desired to determine the steady state phase error 
and the mean time to phase lock for some initial phase offset for the DPLL operating 
with an additive white gaussian noise input. The determination of these two values 
can be made using the same basic test configuration shown in Figure 6-2. This test 
configuration allows the use of identical frequency but phase shifted signals for the 
DPLL's input and reference source. Also, provision is made for the addition of 
gaussian noise to the input signal. 

Referring to Figure 6-2, the source used for both the DPLL reference 
source and the DPLL phase shifted input was a Sulzer temperature controlled crystal 
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oscillator operating at 100 KHz. Since the reference source for the DPLL is re- 

0 

quired to be /’ x f^/ this sets the quiescent frequency of the DPLL at approximately 
390 Hz. To develop an input signal for the DPLL the 100 KHz hardlimited source is 
applied to a presettable 8-bit counter. Then by presetting the counter to various 
values a phase shifted signal at an f rate is obtained. 

The output of the digital phase shift network was then filtered using a 
standard state variable active bandpass filter with center frequency of f Hz and 
Q ~ 100. Utilization of this particular filter design allov/s for easy and precise turn- 
ing of the fi Iter for zero phase shift. The sinusoid output of the bandpass filter was 
then selectably attenuated using Hewlett-Packard Model 355B attenuators to vary 
the signal amplitude. Gaussian noise from a General Radio Type 1381 noise source 
wjs then added to the attenuated phase shifted sinusoid. The power spectrum for 
the Type 1381 noise source is flat over a 25 KHz bandwidth and will therefore look 
to be a white noise source to the DPLL. Noise power was monitored by a HP true-rms 
voltmeter so that accurate signal-to-noise ratios could be set. The corrupted signal 
\\as then hardlimited using a two stage amplifier followed by an LM319 analog 
comparator. Each amplifier stage consisted of an AD518 op-amp configured for 
20. dB of gain at a 40 KHz bandwidth. The output of the LM319 was then used as 
the noise corrupted hard limited input to the DPLL. 

The first test performed on the DPLL involved only the phase detector of 
the loop. For this, the inputs to the divide-by (M xN) and divide-by L counters 
of the loop were removed so that the phase of the DPLL's reference clock could not 
be altered. In this manner, any phase shift between the DPLL's input and reference 




clock that was preset by the phase shift network would be maintained. The proba- 
bility of a count down phase detector output as a function of phase shift and signal- 
to-noise ratio was then determined simply by counting the number of count down out- 
puts and the total number of phase samples taken. The probability of a count down 
output is then the ratio of the former number to the latter, and gives an experimental 
verification of the values of q. as defined by (4-17a). The experimental results are 
plotted in Figure 6-3 and 6-4 and show close agreement with the calculated values 
in all cases. 

The DPLL was then reconfigured for proper operation so that steady state 
phase error test could be performed. For these tests, the phase shift network of 
Figure 6-2 was set to zero so that in the absence of noise the phase error of the 
DPLL would be only the quantization error of the loop The phase state of the DPLL's 
reference clock was then recorded on a Kennedy incremental digital recorder for 
each sample of the input signal. The steady state probabi lity for each reference 
clock state is just the ratio of the number of occurrences of each reference clock 
state to the total number of reference clock samples taken. The variance of the 
reference clock phase is then determined from (4-6). The standard deviation of the 
loop's phase error as determined experimentally has been presented in Chapters IV 
and V along with the theoretical values predicted by the Markov chain model. Each 
of the experimental data points represents at least 300,000 samples of the loop phase 
error and it is obvious that good agreement exists between the experimental and theo- 
retical results. 

The final objective for the hardware tests was verification of the predicted 
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loop transient response. For this, on initial phase offset was applied to the loop 
via the phase shifting network and the number of samples required to the first occur- 
rence of the phase lock state as defined in Chapters IV and V was recorded. This 
was done by detecting the states which defined phase lock as given in Chapters IV 
and V. The detected phase lock state was then used to gate off an event counter 
that was counting to number of phase detector samples taken following the initial 
phase offset. Again, the experimental data has been plotted previously along with 
the theoretical data in Chapters IV and V. Each experimental data point represents 
at least 500 trials and once again close agreement is seen between the theoretical 
and experimental results. 
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CHAPTER VII 


CONCLUSIONS AND RECOMMENDATIONS 

FirsL and second order oil diglLol phase- locked loops have been analyzed for 
both ideal and additive gaussian noise inputs. In addition, a hardware DPLL capable 
of either first or second order operation has been designed and tested for verification 
of the analytical results. For all cases tested, the experimental data showed close 
agreement with the analytical results indicating that the Markov chain model for first 
and second order DPLL's given in Chapters IV and V are valid. 

In Chapter III, ideal inputs were considered for both first and second order 
DPLL“s with the objective of classifying the time response of the loops. For both 
loops it was found that the phase error response was given by a non-linear difference 
equation for which no direct solution was found. However, partial response character- 
istics of the phase error was determined for both first and second order DPLL's when 
the frequency of the input signal is Identical to that of the loop's quiescent frequency. 

Also, expressions for the frequency range over which phase lock can be achieved for 
a fl'st order DPLL and for the frequency range for which a second order DPLL can 
achieve phase lock in minimum time were derived. In both cojses it was found that 
the frequency range was directly dependent upon the number of distinct phase states 
of the reference clock. As would be expected, it wob found that the frequency range 
for which a second order DPLL will achieve phase lock, even with the constraint of 
minimum time to lock, is significantly greater than the frequency range over which a 
first order DPLL will achieve phase lock. 

Specific first and second order DPLL's were also analyzed for stochastic inputs 
by means of a Markov chain model in Chapters IV and V, respectively. From this 
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Markov chain model, the steady state phase error and mean transient response were 
determined. The loop configurations used for the noise analysis were specifically 
chosen to both match the general loop model given in Chapter III and to be realized 
in hardware by standard binary logic families. For both first and second order loops 
it was found that the usual tradeoff between steady state error and transient response 
existed. That is, the steady state error con be decreased only with the cost of a 
longer transient response and the transient response can be decreased only with an 
increase in steady state error. 

For the data presented in Chapters IV and V several specific points are worth 
noting. First, in comparing the transient response of the first and second order DPLL 
it is found that the first order DPLL will achieve phase lock in less expected time than 
the second order DPLL for initial phase offsets less than approximately it/8, Also, for 
a second order DPLL, the steady - t /te phase error degrades rapidly as the signal^to- 
noise ratio decreases below O.U dB. Thus for an application such as an Omega receiver, 
if the initial phase error as the received signal is gated on is expected to be small then 
a first order DPLL will perform In a superior manner over a second order DPLL, However, 
if the signal-to-noise ratio can be expected to be in excess of 0.0 dB and the initial 
phase error is unknown then the second order DPLL will give superior performance. 

Completion of the goals of J'his paper also points to areas in which further research 
would be of use. In particular, for the case of an ideal input, it would be of great 
utility if the determination of the phase error response characteristics were extended to 
include the case of frequency offsets. This extension would then allow the determina- 
tion of time required to lock for an initial frequency offset without the necessity of 
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pe forming a simulaHon. For the case of stochastic inputs, further research on several 
po nts is recommended. First, the concepts of Chapters IV and V should be extended 
to determine pertinent statistics on loop cycle-slippage. Secondly, it would be use- 
ful to develop a comparison of the DPLL's transient response to the standard loop band- 
width used to characterize APLL's. Finally, because of the transient response/steady 
state error tradeoff previously mentioned research into adaptive loops could provide an 
optimum relationship between these two parameters. 
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A. Markov Chains . A Markov chain can be characterized by a system con- 
toining a number of distinguishable states (finite or infinite) for which transition to a 
new state denoted s. from any present state s. depends solely on the present state s.. 
For the application of interest in this paper, assume that the number of system states 
is N, finite. For each state to state transition, define p.. as the conditional proba- 
bility of a transition to state s. given that the present state is s.. Then for defined 

as the probability that the system is initially in state Sj^, the conditional probability 

that the system has traversed through states s., , s._, . . . s. , to the present state 

jl 1^ |n-l 

s.^ is given by. 


P [s. 

r |n 


’jl' ^j2' • • * ^jn-l^ °ilPjr \2^\2' J3 • • • P[n-1, ^ 


For each system state s.,i ~1,2, . , ,,Na vector of transition proba- 


bilities, P. may be written as. 


^i ^*^11 P|2 ■ ■ ‘ PiN^ ' 1 , 2, . . . , N 

Note that the N transition probability vectors must satisfy the requirements. 


Pj. > 0 i, j = 1 , 2, . . . , N 


(A-2) 


(A -3) 


and 


N 


p.. = 1.0 i = 1 , 2, . . . , N 


(A -4) 




and are called stochastic vectors. The N vectors may be arranged in a stochastic 
matrix of size N X N to give the matrix of transition probabilities, [P], 
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P]1 P]2 


P21 P22 


[P] = 


(A-5) 


I^Nl ^^N2 ... 

This matrix then defines all of the state to state transitions for the system and together 
with the initial state distribution a., i = 1, 2, . . N for states s^, ^21 ' ' • 

c(.mpletely defines the Markov chain for states s^, ^ 2 , . . . s^^. 

The matrix of (A-6) defines the probabilities for state to state transitions 

fo‘ a single step. Next consider a transition from state s. to state s. that occurs in 

I ' 

(n) 

exactly n steps . Denote p.. as the probability of observing the system instates. 

at time r + n given that at time r it 'was observed in state s.. Obviously there will 

exist a number of different paths for which the transition from s to s. will occur in 

i I 
(n) 

exactly n steps and p.. ' will be the sum of the probabilities for all of the possible 
paths. For example. 


(A-6) 


k = 1 


Pik Pkj 


(A-7) 


It can be shown by induction that p.. is given by the recursive relationship, 


(n) ' (n-1) 

Pki 


(A-8) 




If the value are arranged as elements of a matrix denoted [P]*^ then It is 

M 

obvious from (A-8) that standard matrix multiplication is applicable so that 

[P]" = [P] [P]""^ (A-9) 

Recall that a. was defined as the probability of observing the system in 
state s. at time zero. Then the unconditional probability of observing the system in 
state s. after n steps is 
N 


a':") == I 


(n) 

a. p.'. 
i = 1 ' 


(A-10) 


(n) 

Further the distribution of a ; ' will tend to be independent of the initial distribution 
' 1 

a. , i = 1,2, . . . , N if p!?^ is independent of i which will be the case if [ P 
converges to a matrix of identical rows. If [ P does converge, then as n -► oo the 
distribution a^*^^ is the steady state distribution tor the system. 

It is now necessary to show that such a steady state distribution exists. 
First, define the states of an aperiodic Markov chain as being transient if the proba- 


bility of reoccurrence of that state is less than one for infinite time so that 

Z pW < a, 
n= 1 " 


(A-11) 


and define the state as being null if the probability of reoccurrence of that state is 
one but the mean time to reoccurrence is infinite. It has been proven by Feller [8], 
Chapter XV. 6 that for an aperiodic Markov chain that all states are either transient 
or null or all states are ergodic. If all states are ergodic then 


lim 

n -► 00 


» - 


kk 


= U. > 0 


(A-12) 


where 1”^® mean reoccurrence time of state Sj^ and Uj^, k = 1,2, . . . N 
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is the unique steady state distribution for the system. 


A Markov chain can also be specified for which each state has one or 
more possible vectors of transition probabilities. For each state s, , the possible 
vectors are called alternative vectors and transition from the present state to some 
new state is governed by one and only one of the alternative vectors associated 
with the present state. In this case, for each state s, there exists k. alternative 
vectors. 


k = (k k . . . k ) i = 1, 2, . . . N 
Pi Pil Pi2 PiN k = 1,2, . . . K. 


(A-13) 


where 'p is the probability that the system will make a transition to state s. given 

'i I 

that the present state is s. and k alternative for s. is used. As before, each of the 
alternative vectors of transition probabilities are stochastic vectoR and must satisfy 
the conditions , 

k > 0 i, i = 1, 2, . . ., N and k = 1, 2, . . . K. (A-14) 

Pii “ 


= 1.0 = 1, 2, . . . , N and k = 1, 2, 


.,K. (A-15) 

I 


The vectors may be combined to give the K x N stochastic matrix, P , 



B. Program Listings . The following three fortran programs were written to 


solve for the DPLL's performance characteristics as described by the Markov model 
given in Chapters IV and V. 

1 . PBSTGEN . Program PBSTGEN is used to establish the systems of 
equations associated with ( 4-1 )/(5-2 ) and ( 4-9)/ (5 -10) . For these systems of 
equations it is seen that only two elements of each row are non-zero. Therefore, 
instead of storing the entire coefficient matrix for these systems of equations, PBSTGEN 
generates two vectors whose lengths are twice the dimension of the coefficient matri- 
ces for each system. The first vector, PNTRi for (4-1 )/(5-2) and T1 for (4-9)/ 
(5-10), contains the alternative loop state number for each equation (i.e. the pos- 
sible state to state transitions for each equation) while the second vector, PNTR2 
for (4-1 )/(5-2) and T2 for (4-9)/(5-10) , contain a pointer for the probability co- 
efficient associated with each state to state transition. The actual value of these 
probability coefficients will be calculated in the following programs. 

2. PBDPI-2 . Program PBDPL2 uses the vectors PNTRI and PNTR2 
generated in PBSTGEN to solve for the steady state phase error of the DPLL. Solu- 
tion is performed by means of Jacobi's iterative technique until the steady state prob- 

- 8 

ability for all alternative loop states converges a difference of less than 10 for 
successive iterations. If convergence is not achieved, a maximum of 10,000 itera- 
tions are performed. The iterative solution is performed for signal-to-noise ratios 
from -40. dB to 20. dB in increments of 5 . dB and the output consists of the steady 
state probability of each reference clock state and the standard deviation of the phase 
error. 
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3. PBDPL2T . Program PBDPL2T uses fhe vectors T1 and T2 generated by 


PBSTGENto solve for the mean time to phase lock for an initial phase offset. Solu- 
tion is performed by Jacobi's iterative technique until the result converges to a dif- 

-5 

ference of less than 10 for successive iterations. Again, a maximum of 10,000 
iterations will be performed. The solution is performed for signal-to-noise ratios 
from -40. dB to 20. dB in increments of 5. dB. The output consists of the mean time 
to lock for initial phase offsets from -tt to ir where the value of the divlde-by L 
counter is initially zero. 
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C. Binary Phase-Locked Loop Design. The following sections discuss in 


detail the design of the hardware DPLL used to validate the Markov chain model. 

The block diagram for the DPLL may be found in Figure 6-1 . 

a. Binary Phase Detector . Referring to the schematic of Figure C-1 , the 

binary phase detector operates by sampling the binary input signal f and producing 

c 

a complemented pulse output on either the count-up or count-down line. This func- 
tion is produced by three D-type flip-flops as follows. The OPEN SW signal is 

applied to the clock input of ff B3 while the binary signal f is applied to the data 

0 

input of the same ff. Thus on a positive transition of OPEN SW the value of f^ is 
latched, giving the sampled value IN at B3's Q output. The two ff's of A6 are 
initially set to the ONE state so that when SET lAT latches the input values to 
the ff's (IN and IN), one ff goes to the ZERO state while the other remains in the 
ONE state. The SET LAT signal is followed by the CLR LAT signal which sets both 
ff's to the ONE state. Thus the count-up or count-down are produced in a mutually 
exclusive manner when one of the ff's of A6 toggles HIGH LOW HIGH while the 
other remains HIGH. 

b. DIVIDE-by L Saturating Counter. The saturating up-down divide-by 
L counter is shown in Figure C-2. The counting function is performed by two ssrlnUy 
connected 74193 4-bit synchronous up-down counters. The output states of the 
counter is detected by a logic network to produce the INHIB UP and INHIB DNI 
signals that will inhibit the UP and DN clock signals respectively. For example, 

if switches SWl through SW6 are closed, then for a counter state of 0000 0001 the 
INHIB UP signal will be TRUE, thus inhibiting the UP clock signal and saturating 
the courjter at that value. Note, a DN clock signal will still count the counter 
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to the 0000 0000 state. Similarly, for all switches closed and the counter In state 


1111 nil the INHIB DN signal will be TRUE and the counter Is saturated at that 
state. If switch SWl is opened while all other switches remain closed, then the 
counter saturation states will be 0000 0011 and 1111 1101 (±3). Continuing in this 
manner, the counter saturation states may be selected to be ±2' - 1 , i = 1 , 2 , . . . 7 
by opening switches SWl through SW ( i- 1 ) whi le all other switches remain closed . 

c. Divide-by (M x N) Counter . The schematic for the divide by 

(M X N) counter along with the buffer register and adder are shown in Figure C-3. 
The counter consists of three serially connected 74193 4-bit up-down counters. 
Following an UP or DN clock signal the value of the divide-by (M x N) counter is 
loaded into the 12-bit buffer register composed of two 74174 hex D-type flip-flops 
by the LOAD2 signal . The values of the 12-bit buffer and the divide-by L satura- 
ting counter compose the inputs to the adder circuit consisting of three 7483 4-bit 
binary adders. The output of the 12-bit adder is applied to the preset terminals of 
the divide-by (M x N) counter. If the type select switch is set for first-order opera- 
tion, the adder output will not be loaded into the divide-by (M x N) counter. If 
a second-order loop is selected, the LOADl signal will be applied to the load in- 
puts of the 74193's thus presetting the divide-by (M x N) counter to the value of 
the 12-blt adder output. 

d. Variable Phase Reference Clock . Referring to Figure C-4, six 

selectable contiguous bits from the divide-by (M x N) counter are applied to one 

input side of a binary magnitude comparator formed by two 7485 4-bit magnitude 

comparators. The other input to the comparator is obtained from the lower N bits 

N 

of the two series-connected 74193 counters that are being clocked at a rate of 2 

-no- 






time the input frequency. Thus the EQUALS output pin 6 of A3'( , of the magnitude 

comparator will be a ONE when the two inputs are equal and this will occur at a 

rate equal the input frequency f . The circuit formed by A33 and A34 prevents the 

0 

sample output from occurring at a rate greater than f . For example, if pin 8 of A34 
is HIGH then when the magnitude comparator detects the equal condition, pin 5 of 
A34 is set HIGH causing the phase detector to sample the input signal. The SAMPLE 
signal is reset LOW by the LOAD1 signal at the end of the loop phase update. How- 
ever, since the (M + N)-bit counter may have counted up during the last phase up- 
date and the counters of the reference clock have also counted up, it is possible for 

N *f 

the magnitude comparator to detect equal states occurring at a rate of 2 c. 

Thus it is necessary for the LOADl signal to reset pin 8 of A33 to a ZERO so that 
another sample cannot occur until the reference clock counter loads a ZERO into 
that flip-flop. 

e. Control Logic . A timing diagram for the necessary control v/ave- 
forms and the logic schematic used in their generation are given in Figures C-5 and 
C-6, respectively. Use of the control waveforms have been described in the pre- 
vious sections. 
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Figure C-5. Control Waveforms. 
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D. Digifal Low-Pass Filter. Consider a divide-by M binary up/down counter 


whose input x(t) = ± 1 is clocked at a rate f = 1/T. The contents of the counter 

c 

will be the sum of the previous M- 1 inputs plus the present input, so that the value 

of the counter y(t) may be expressed by the discrete equation 
. M- 1 

y(kT) = i Z x[(k-i)T] (D-1) 

i = o 

up to the point of counter overflow or underflow. Taking the Fourier transform of 
( D-1 ) gives 


1 


M- 1 


Y(M=j^Z 


I = O 


which gives the transfer function for the counter as, 

. AA — 1 • • 

u/* \ ' V ..“ 1 '“^ 

H(p) = I 

I = o 


i 1 - e 

M 


-juTM 


1 - e 
sin / uTM' 


-juT 


I \ 2 I -juT(M-l) 

M . ,oT . ^ 2 

sin 


(D-2) 


(D-3a) 


( D-3b) 


(D-3c) 


Substituting for T = 2ir/u in (D-3c) gives the mcsgriitude-squared function for the 

c 

counter as , 



2 r 

H,. . 
(lu) 

II 

1 — 


1 sin (ttM (j/(jc) 
M sin (it (j/uc) 


(D-4) 


which has been plotted in Figure D-1 for M = 2,4, and 8. As can be seen from 
Figure D-1 , the binary up/down counter acts as a digital low-pass filter whose band- 
width decreases as M increases . 


- 118 - 



